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Abstract. We give a geometric proof of spectral stability of travelling kink wave solutions to 
the sine-Gordon equation. For a travelling kink wave solution of speed c ^ ±1, the wave is 
spectrally stable. The proof uses the Maslov index as a means for determining the lack of real 
eigenvalues. Ricatti equations and further geometric considerations are also used in establishing 
stability. 
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^ ' 1. Introduction 

The sine-Gordon equation: 

u tt = u xx + smu (1) 

has applications in many areas of physics and mathematics. It can be used to model magnetic 
^ ■ flux propagation in long Josephson junctions: two ideal superconductors separated by a thin 
insulating layer [SCR76J, |DDvGV03] . It can be thought of a model for mechanical vibrations 
of the so-called 'ribbon pendulum' - the continuum limit of a line of pendula each coupled to 
their nearest neighbor via Hooke's law [BM07] . In biology, it has found applications in modeling 
the transcription and denaturation in DNA molecules [Sal91j . Further, it can be used to model 
propagation of a crystal dislocation, Bloch wall motion of magnetic crystals, propagation of 
"splay-wave" along a lipid membrane, and pseudo-spherical surfaces to name a few others (see 
|SCM73] and the references therein). 

In this paper, we consider solutions of the form v(x + ct,t) where c is the (positive) speed of 
the traveling wave. Making the change of variable z = x + ct and substituting into (JT|) gives 

c 2 v zz + 2cv zt + v tt = v zz + sin v. (2) 

A travelling wave solution will be a t independent solution v(z) to fl2]). Thus it solves the 
(nonlinear) pendulum equation: 

(c 2 - l)v zz = sin v. (3) 

A kink wave solution is a travelling wave solution to (J2]) corresponding to a heteroclinic orbit in 
the phase plane of d3J). 

Stability of a singularly perturbed kink wave solution was shown in [DDvGV03j . In [BM07J, 
the Cauchy problem for the sine-Gordon equation in laboratory coordinates was studied using 
inverse scattering techniques. In this paper we take a more geometric approach determining 
the spectral stability of kink-wave solutions to equation ([1]) via geometric considerations and 
elementary methods of ODE theory. 
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1.1. The Maslov index. Note: In this paper we will use the description of the Maslov index 
as in [ESS]. 

Let the matrix J = ^ denote the standard symplectic structure on IR 2 . A line £ passing 

through the origin is considered Lagrangian in the sense that for any V\ , v<i G £ the inner product 
of v% with Jv2, < Vi, Ji>2 >= 0. Let £{t) be a curve of lines in R 2 . If £(t) can be written as 

fx(t)\ 

the span of the 2x1 matrix I y)Aj> we wm can the functions x[t) and y(t) a frame for £{t). 

Alternatively we can view £(t) as a curve in RP 1 w S , and if ( J * s a f rame °f ^? then 

£{t) = [x(t) : y(t)]. Now let a = [a x : a 2 ] be a fixed line in R 2 . Suppose that £{t) : [t ,* n ] -> MP 1 , 
and £{t) — a at ti,tz, ■ ■ ■ t n -i with t < t\ < t 2 < ■ ■ ■ < £ n _i < t n - Suppose further that we have 
£(ti) 7^ (0,0) for all U. We define the crossing form, T(£{t),a.,ti) of £{t), with respect to a, at U 
as: 

a, u) = x{ti)y(ti) - y(ti)x(ti)- (4) 
For a curve £{t) in MP 1 as above, we define the Maslov index, fi(£(t) , a) , as: 

M^),a) = 5>gn(r(^),a). (5) 

The Maslov index defined in this way is a signed count of the number of times that £(t), viewed 
as a curve in S" 1 crosses the point a. 

2. Kink Waves 

Travelling wave solutions to the sine-Gordon equation for which the quantity c 2 — 1 < are 
called subluminal waves. When c 2 — 1 > they are called superluminal waves. We have the 
following theorem: 

Theorem 1. Kink wave solutions to equation (TJP u tt = u xx + sinu, are spectrally stable if c 2 ^ 1. 

2.1. Subluminal kink waves. Until otherwise specified, the quantity c 2 — 1 < 0. We will focus 
primarily on the orbit that satisfies the following boundary conditions, though the analysis that 
follows will apply to all subluminal kink waves. Let v(z) be a solution to (j3J) satisfying: 

lim v(z) = —tv and lim v(z) = n (6) 

Linearizing equation §2§ about the kink wave solution v gives: 

(c 2 - l)(f zz + 2ap zt + iptt = (cosv)ip. (7) 
Setting <p t = if; we can rewrite (j7j) as: 

'<P\ = $ 

J t cos vip - (c 2 - l)tp zz - 2cif) 

In order to establish spectral stability, we need to consider the eigenvalue problem of equation 
(jHl). Letting A be the eigenvalue parameter, we obtain: 

^ = A 

(cosv)ip - (c 2 - l)(f zz - 2cip z Uf) 



,2 ^ o-./. (8) 



2 a. o«/. - A r, (9) 
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We will consider (twice differentiable) L? perturbations, (p. This leads to the following eigenvalue 
condition: A function ip is an eigenfunction of the linearized operator with eigenvalue A if 

(c 2 - l)ip zz + 2c\(p z + (A 2 - cos v)<p = 0, (10) 

is satisfied, together with the condition that 

lim <p(z) = 

z— >±oo 

The idea now is to reformulate the above eigenvalue condition in an equivalent, geometric way. 
Setting w = (wi, w^), with w\ = cp, wi = tp z and ' = we have: 

w ' = (!? J = (cos^-a 2 -2cA J L? J =:A(\,z)w. (11) 



Also we set 

A(\):= lim A(X,z)= (J_ x * _LV (12) 

We remark that A(X) has an unstable and a stable subspace which we will denote by £ u and 
£ s respectively. 

Lemma 2 (Geometric version of the eigenvalue condition). Given that (^M solves (TT\) 



is equivalent to 



lim pi = (13) 

2— >-±0O 



lim ( Wl J C, and lim ( Wl ] -> f (14) 



Proof. We compactify the extended phase plane of (TTTT) by introducing a new variable z = z{r) = 
tan(^r), where r G (—1, 1). Now the extended system of ODE's becomes the (autonomous) 

w' = A(t,X)w , v 

7T Z 

Next we note that if we view the space 1R 2 x (— 1, 1) as a plane bundle over the segment [—1,1], 
the linearity of ( TTTj) means that linear subspaces of the fibres are preserved, so the flow defined 
by flUD onl 2 x (-1,1) induces a well defined flow on the cylinder S 1 x (—1, 1). Moreover, this 
flow can be continuously extended to a flow on [—1,1]. Finally we note that on this cylinder, 
the points £ u and £ s are fixed points of the induced flow, with a one dimensional unstable and 
stable manifold respectively. The uniqueness of the stable and unstable manifolds concludes the 
proof of the lemma. □ 

We have redefined the original eigenvalue condition as the existence of a heteroclinic orbit 
joining and £ s . 

Now, as c 2 — 1 < 0, the matrix 

A(\,z)->A(\)=(_Zx» ieA 
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has eigenvalues j u and 7^ corresponding to the unstable and stable subspaces. 



<y - ~ CA - ^ ~ ^ ~ B and 7 - Z^WgEEES fl6) 
<r — 1 <r — 1 



and (c 2 - 1)7 2 s + 2cA7„ iS + (A 2 + 1) = 0. Observe: 



' -c\+-Jx 2 -(c 2 -l)\ /-cA- A /A 2 -(c 2 -l)' 

f 4 =< I a^Ti > , and C =< —1 I > . (17) 



where < ( > denotes the linear space spanned by the vector y^j . For each fixed A, let £(z) 

be the set of lines in R 2 that tend to the unstable subspace of A at —00, under the flow of ( TTT1) . 
that is: 



£( z ) = <j < >\(Z 1 ) solves J10J and ->■ as z -> -00 }> (18) 



The lines £(2) are Lagrangian, and in fact the representation L^M ^ s a f rame f° r the line so 

we can define the crossing form relative to the subspace £ s , F(£(z), £ s , A, z) as T(£(z),^ s , A, z) = 
Wiu^ — W2w[. Substituting as in equation ffTTj) gives 



^/^/ , , c A . cosw — A 2 9 2cA 9 . . 

r (£(z , A, 2 = - w\ - -^—-w x w 2 - w\ 19 
cr — 1 c z — 1 



Lemma 3. T/ie crossing form of£(z), relative to the stable subspace at infinity, T(£(z), £ s , A, z) 
defined above is independent of X. 
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Proof. To evaluate the crossing form on £ s , we use the fact that on £ s we have that W2 

A 2 + l (A 2 + l) rp, • 

/ Wi = v , a ^ i. lhis gives 

-c\+y/\2-(c2-l) 7a C 2 -l) & 



cost; -A 2 2cA (A 2 + 1) / (A 2 + 1) 



/ 1 2 


(c 2 - 


1) 






w ? 




(c 2 - 


1) 










(c 2 - 


1) 




9 


, w l 


(c 2 








u\ 




(c 2 - 


1) 






w ? 




(c 2 - 


1) 



c 2 — 1 c 2 — 1 7s (c 2 — 1) V 7 S (c 2 — 1 



(cost; - A 2 ) 7 2 (c 2 - 1) - 2cA(A 2 + l) 7s - (A 2 + 1) : 



7l(c 2 - 1) 

cost^c 2 - 1) - (A 2 ) 7 2 (c 2 - 1) - 2cA(A 2 + l) ls - (A 2 + l) 2 

7 2 (c 2 - 1) 

cost;7 2 (c 2 - 1) - (A 2 ) 7 2 (c 2 - 1) - 2cA(A 2 + l) 7a - (A 2 + l) 2 



7 2 (c 2 - 1) 
-7, 2 (c 2 -l) + 7, 2 (^ 2 -l) " 

,,2^2 



7 2 (c 2 - 1) 

(cost; + l) 7a V - 1) - (A 2 + l)((c 2 - 1) 7 2 + 2cA 7s + (A 2 + 1)) 

7 2 (c 2 -l) 



(cost; + l)7 2 (c 2 - 1) 



7 2 (c 2 - 1) 

because (c 2 — 1) 7 2 — 2cA7 s + (A 2 + 1) = 0, so we have that 

(cost; + l)ujf 



Which shows that T is independant of A. 



(c 2 - 1} 



(20) 
□ 



Evaluating the limits and using lemma |3] we get that the Maslov index of the curve of La- 
grangian subspaces £(z), with respect to the stable subspace at infinty is: 



z 

= — # of crossings that occur. 



(21) 
(22) 



We next observe that only regular crossings occur. A crossing is called regular if &(z) is non- 
zero. Equation (T2"Uj) is clearly non-zero, except at z — ±oo. 

At a regular crossing we have that £(z,X) = = £ S- ^ ne f unc tion f(z,X) = is well 

defined at a regular crossing of the unstable manifold £(z) with £ s . Moreover at a regular crossing 
we have that 

„ .x A 2 + l 
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and that 

a /M = WAi ^ (24) 

Thus we can use the implicit function theorem to write the location of the crossing in the z 
variable as a function of the parameter A, z = z(X), further because equation (I20I) is independant 
of A, we can do this for all A. 

When A> 1, the system in ([TT]) tends to 



A 2 -2cA 



^2/ V^23T ^231/ \«W 

—A 



= A ( X » 1) ( 25 ) 



which has eigenvalues y± = Which are real and as c — 1 < 0,one is positive and one is 
negative. For large enough A then, the number of crossings of £{z) relative to £ s is the same 
as the number of crossings as for the solution of the constant coefficient equation (1251) tending 
towards the unstable subspace of the matrix A(X ^> 1). Thus for large enough A the number of 
crossings of the unstable manifold £(z) and £ s is equal to zero. 

Now fix a A and a location of a crossing Zi(X) say. Because the number of crossings tends to 
zero as A increases, and because of the implicit function theorem assertion above, we have that for 
some finite Ai, lim Zi(A) = oo. But this is exactly the geometric reformulation of the eigenvalue 

A— >Ai 

condition. Moreover, the implicit function theorem implies that we have a unique value, ^i(A), 
for this crossing for each A, and the locations of different crossings (say .22(A)) cannot intersect. 
Thus as A increases, the number of crossings will decrease monotonically, as will the Maslov 
index. This last paragraph is summed up in the following corollary. 

Corollary 4. If N\ is the number of times £{z) crosses £ s when X = Xi, as z ranges from — oo 
to oo and N 2 is the number of times £(z) crosses £ s when X = X 2 , then | A^! — iV 2 1 = the number 
of eigenvalues of the linearized operator (TJj) in (Ai, A2). 

Lastly, we show that p,(£(z), £ s , 0, z) = 0. This shows that there are no crossings when A = 
and so there are no crossings for A G (0, 00). To see this observe that when A = 0, equation (fTTj) 
reduces to the equation of variations of the pendulum equation ([3]). Thus the unstable manifold 
of equation (fTTj) is just the tangent line to the heteroclinic orbit, joining (— n, 0) to (71", 0) in the 
phase plane of (jllD . This orbit in the phase plane is given by: 



2^±i (26) 



c 2 



and the slope of the tangent line is given by 

dv z — sin v 



(27) 



dv ^(c 2 - l)(2cosu + 2)' 

In order for a crossing to occur, the tangent line must be parallel to £ s . This means that 

— sinv —1 

V / (c 2 -l)(2cosw + 2) ~ ^-(c 2 - 1) 

But this only happens when v = ±7r. But v = ±7r is a critical point in the phase plane. 
That is, there is no z G (—00,00) where this can happen, so there can be no crossings, so 
n(£(z),£ s , 0, z) = 0. This shows that there are no crossings for A G (0, 00). 



(28) 
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Lastly we note that the crossing form argument works in the negative A direction and that the 
asymptotic behaviour of the system in (fTTj) for A — 1 is the same as that for A ^> 1. Thus we 
have no real eigenvalues A ^ 0. 

2.2. Complex Eigenvalues. We now return our attention to equation ffTUl) to show that there 
are no complex eigenvalues A = p + iq where p > . Without loss of generality we can assume 
that q > so A = re 10 where r > and 9 £ (0, |). Rewriting ffTOj) as: 

„ , 2cA , (A 2 - cos(t;( g))) 

^ — i" + 2 — r 

— 1 CT — 1 



w" + — V + * \ K,,, w = 0, (29) 



we make the substitution 
and re-write (1291 as 



^ = we te z ) (30) 



c 2 -l (c 2 -l) 

We wish to consider now complex valued solutions to (13"T|) . We let 77 = ^ and consider the 

induced flow on a chart of CP 1 , where 1/) 7^ 0: 

, _ ^> - y>y _ cos v - A 2 c 2 A 2 2 

77 " ^ " c 2 -l + (c 2 - l) 2 V • (d j 

If we set rj = a + and consider the imaginary part of the vector field of ( |32l) on the real axis, 
that is those points where rj = a, we have 

^ l/3=o = + (c 2 _ 1)2 = (c 2 _ 1)2 > ifp,?>0. (33) 

thus the flow of fl32|) on the real axis of CP 1 is always pointing in the positive imaginary direction. 

By considering ^ from equation (fTTj) . we now interpret the eigenvalue conditions, for the 
eigenvalues that we are interested in. 

Lemma 5. An eigenvalue X is a value of A where there exists a heteroclinic orbit of the flow 
induced by on (a chart of) CP 1 , i. e. A an eigenvalue means that there is an orbit from 

lim — = 7„ to 7 S = lim — . Under the transformation given in we have that 

W2 cA 

— + -2—1 = 

Wi c z — 1 
and so 7„ goes to rj u and 7^ goes to r) s where 

V A 2 - (c 2 - i; 
c 2 - 1 
y/A ^T^ l) 
c 2 - 1 

We claim that for A = p + iq with p and q both positive, that r/ u has a positive imaginary part 
and i] s has a negative imaginary part. Thus as /3'L_ > 0, from above, there can be no orbit 
joining the two critical points, and so no eigenvalues with positve real part. Thus all kink waves 
satisfying c 2 — 1 < are stable. 



V=- + —> (34) 



Vu = — — ~ 2 , and (35) 
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To see that r] u has a positive imaginary part, write A = re ld and r\ u = R u e t0u . Since we are 
assuming that A has positive imaginary part 6 G (0, |). By equation we have that 

/ r 2 sin(2#) \ 

9 U = arctan — — . 37 

\r z cos(2#) + (1 — cr) J 

Since 9 U will be the same for all points on the line re t9u , and since c 2 — 1 < 0, we can assume 
without loss of generality that r 2 = (1 — c 2 ). Thus substituting into equation ( 1371) we have: 



/ sin(2fl) \ / 2sin6>cos6> , 

U u = arctan — = arctan — — — 5— = 0. (38) 

1 cos(2#) + 1 J \cos 2 (6) - sm 2 (9) + 1 ' 



The exact same calculation shows that rj s has a negative imaginary part. 

We have shown that there can be no homoclinic orbit on the chart of CP 1 parametrized by 

— . To see that on there can be no such orbit on the other chart, we have that off the origin, 

nil 101 

the charts are transformed into each other via rj — > - for £ G C\0. Thus r\ u — V — which will 
have a negative imaginary part. Likewise — will have a positive imaginary part. If we write 
( = = d + iq and consider the flow on this chart induced by ([10]) 

_ (t//) 2 - W _ , cosv - A 2 c 2 A 2 
C 1 - ( C 2_! + j^YyK ■ (39) 



Now letting C = a + ir and considering only the imaginary part of the flow on C given by ([39 
restricted to where r = we have that: 



r L=n= —2 -j—- 7 = - 772 U2 <0 ifp,g>0. (40) 



-2pqc 2 cr 2 2pqa 2 2pqa 2 

(c 2 - 1) 2 ~ ^f^l = ~{c 2 - \y 

Thus the flow of (|32|) on the real axis of (this chart of) CP 1 is always pointing in the negative 
imaginary direction, and so there can be no heteroclinic orbit connecting — to — , and hence 
there are no eigenvalues A with positive real part to equation (jUJ). 

The same argument can be run in for A = p + iq, with p and q both negative to show that 
the flow on the real axis is pointing in the wrong direction to allow for an orbit joining 77" to 
r] s as well. Just as in the previous case, one makes a Liouville transformation and then tracks 
the location of the stable and unstable subspaces on both charts of CP 1 . Thus we conclude that 
there are no eigenvalues off the imaginary axis. 

2.3. Superluminal kink-waves. We now return our attention to the superluminal eigenvalue 
problem equation <M): 



2 o™, - A S (41) 



^ = A ( V 

(cos v)(p- (c 2 -l)(f zz - 2cip z \ip 

We are interested in the values of A G C for which there are solutions to (14TT) which decay to 
zero as z — > ±00. We first investigate the limiting cases lim^-too, which becomes the constant 
coefficient ODE 

(1 - c 2 )tp" - 2c\ip - (A 2 + !)<p = (42) 
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c 2\ 



which can easily be solved for any value of A G C. The characteristic exponents are given by: 

cA ± y/A 2 + ( 1 - 
1 -c 2 

By taking limits asc-f 1 + , we obtain that the signs of the real parts of r\p are equal. In fact 
they will be the opposite of the sign of the real part of A. We have: 



ri, 2 = w - i (43) 



sgn ( dt( lim ri) ) = — sgn (5ft(2A)) and, sgn ( dt( lim r 2 ) ) = — sgn (dt(\ + — )) 



Moreover, by differentiating each in the c variable, we obtain that the sign of the derivative 

>1,2 

dc 



of the characteristic exponents is the same as the sign of the real part of A. Thus we can 



conclude that in the case of c 2 — 1 > 0, we have that for the asymptotic cases, the eigenvalues 
of A(X) have the same sign. This means that as z — >■ — oo there can be no unstable orbit of 
the construction in lemma [2] and therefore there can be no heteroclinic connection to the stable 
orbit at oo. Thus we conclude that there is no point spectrum when c 2 — 1 > is positive. This 
concludes the proof of theorem [1] 
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